###Exerccio 4.2: Cervejas###
cervejas <-c("chope","lata","garrafa","chope","garrafa", "garrafa","lata","lata","nenhuma","lata","garrafa","garrafa", "garrafa","lata","lata","lata","garrafa","lata","chope","nenhuma", "garrafa","garrafa","garrafa","chope","garrafa","garrafa","chope","garrafa","lata","lata")
#Grfico de barras
barplot(table(cervejas))
#Grfico dotplot
dotchart(table(cervejas))

###Exerccio 4.3###

#1.Histograma do dao dos fustes dos caixeitais
caixeta=read.table("caixeta.csv",header=T, sep=",")
str(caixeta)
summary(caixeta)
class(caixeta)
hist(caixeta$cap/pi)

#2.Histogramas da altura das rvores para diferentes caixatais ("local")
par(mfrow=c(3,1))
hist(caixeta$h[caixeta$local=="chauas"])
hist(caixeta$h[caixeta$local=="jureia"])
hist(caixeta$h[caixeta$local=="retiro"])

#3.H diferenas entre as estruturas (distribuio de tamanhos) dos caixatais?
#As medias e distribuies diferem no apenas entre os tamanhos, mas tambm entre a frequencia de ocorrencia em cada local.

### Exerccio 4.4 ###
e.grandis <- read.table("egrandis.csv",header=T,sep=";")
summary(e.grandis)
class(e.grandis)
head(e.grandis)

#1. Utilize o grfico boxplot para analisar o DAP de rvores de E. grandis em funo das variveis regio (regiao) e rotao (rotacao). 

par(mfrow=c(2,1))
boxplot(dap~regiao,data=e.grandis)
boxplot(dap~rotacao,data=e.grandis)

#2.Avalie a normalidade da altura do conjunto total de rvores com um grfico quantil-quantil contra a distribuio normal. 
par(mfrow=c(1,1))
qqnorm(e.grandis$ht)
qqline(e.grandis$ht)
##A varivel altura nao segue uma distribuio Normal.

### Exerccio 4.5 ###
#1.Analise a relao dap-altura ('dap' e 'h') em funo do caixetal (local) com a funo plot, mas somente para as rvores de caixeta (Tabebuia cassinoides). 
caixeta<-read.table("caixeta.csv",header=T, sep=",")
summary(caixeta)
class(caixeta)
head(caixeta)
caixeta1<-caixeta
caixeta1
caixeta1$dap<-caixeta1$cap/pi
caixeta1
head(caixeta1)
nome.spp<-unique(caixeta1$especie)
nome.spp
caixeta1[caixeta1$especie==nome.spp[3],]
par(mfrow=c(3,1))
plot(caixeta1$h~caixeta1$dap,data=caixeta1,subset=local=="chauas")
plot(caixeta1$h~caixeta1$dap,data=caixeta1,subset=local=="jureia")
plot(caixeta1$h~caixeta1$dap,data=caixeta1,subset=local=="retiro")

#2.Verifique a linearidade com a funo scatter.smooth
par(mfrow=c(1,1))
?scatter.smooth
scatter.smooth(caixeta1$h~caixeta1$dap,span=1/2)

#3. pacote lattice
library(lattice)
xyplot(caixeta1$h~caixeta1$dap|local, data= caixeta)

